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ABSTRACT 


The subject of this paper is an investigation of the 
significant factors that affect the performance of a sub- 
cavitating pitot type intake for a waterjet propulsion 
system for hydrofoil craft. A computer model for design 
of an inlet nacelle is developed which accepts values of 
two parameters and several system constants from a control 
program and designs a nacelle to meet those requirements 
subject to limitations imposed by cavitation and geometry. 
The parameters are inlet velocity ratio and inlet diameter 
to maximum diameter ratio. A body of revolution config- 
uration with profile shape based on certain desirable 
pressure distributions is assumed. Extensive use is made 
of empirical data from the available unclassified literature. 

The design conditions are those existing at maximum 
speed, called cruise speed, and the minimum foil-borne 
Speed, called take-off speed. 

The subroutine will accept up to three additional speeds 
at which a performance prediction will be made after design 
is completed from th first two. 

Pressure recovery characteristics and cavitation limits 
for some typical designs are displayed as graphs. 

The subroutine listing and flow diagram is included in 
Appendix C. In addition, an alternate design method which 
accepts only the inlet velocity ratio as a parameter is 
included in Appendix D. 
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Chapter 1. 
PROBLEM STATEMENT 

Ram (pitot) type inlets for waterjet propulsion systems 
in hydrofoil craft are responsible for a significant part 
of both the external drag and internal flow losses and 
have a strong influence on the propulsive efficiency of 
the system. 

The objective of this investigation is to identify and 
quantify the significant characteristics of a "well de- 
signed" inlet in a manner that will allow their use at the 
preliminary design stage to assist in rapidly determining 
the feasibility and relative merit of various propulsion 
system configurations in a craft of a given size with 
given speed and range requirements. 

The term, "well designed", as used with reference to 
the inlet and nacelle will mean that inlet configuration 
which best suits the propulsion system as a whole when 
optimized with respect to some chosen criterion. 

À computer subroutine is to be developed for use with 
an optimization program concurrently being developed by 
R. P. Gill (ref. 1) which calls this subroutine and others 
written by Gill and R. C. Percival (ref. 2). The optimizat- 
ion criterion chosen for this particular program is least 
propulsion system weight including fuel, but the inlet 
subroutine, called SUBROUTINE NACEL, does not depend on 
that particular criterion and could be used for others. 


Engineering approximations and empirical fits are 
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used wherever they are suitable and consistant with the 
purpose of the study. 

There are at least three broad categories of nacelles 
which could be used for ram inlets: subcavitating, super- 
cavitating, and base-vented. Only subcavitating nacelles 
are considered in this study. Persons who are interested 
in supercavitating and base-vented nacelles may wish to 
read references (3) and (4) for further discussion on 
those types. 

Subsonic aerodynamics provides much of the data 
utilized in this study, but hydrodynamic design has the 
additional complication of cavitation avoidance. For 
the purpose of this paper cavitation is assumed to start 
when the local static pressure decreases to the vapor 
pressure of sea water. Incipient cavitation is imposed 
as an absolute limit on the decrease in static pressure 
at the design conditions even though a certain amount of 


cavitation may be allowable in a real life situation. 
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Chapter 2. 
DESIGN CONDITIONS 

There are three possible flow conditions at the inlet 
and Bernoulli's equation indicates the changes in local 
Static pressure coefficient associated with each of them. 

When Vi/Vo<1, the stream lines are diverging (Fig.1) 
and since there are no losses (to the inlet plane, at 
least), local velocities inside are lower than free stream 
velocity causing the pressure coefficient to be positive 
and no cavitation danger exists. Local velocities on the 
exterior of the forebody are higher than in the free stream 
and the pressure coefficients are negative (Fig.2), so 
Cavitation danger exists and increases as Vi/Vo decreases. 
Thus, for any particular nacelle and free stream velocity 
combination, there is a minimum velocity ratio below which 
Cavitation occurs in the external flow. 

When Vi/Vo>1, the stream lines are converging 
(Fig. 1). Local velocities in the external flow will 
Still be somewhat greater than in the free stream but 
lower than in the diverging flow case, so if the external 
flow was cavitation free when Vi/Vo <1, it will also be 
cavitation free when Vi/Vo>1 at the same Vo. Local 
velocities in the internal flow are now greater than the 
free stream with a corresponding decrease in static pressure 
coefficient, and cavitation inside the lip is a definite 
possibility. There is therefore a maximum velocity ratio 


for each nacelle and Vo above which cavitation will occur 
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Figure 1. Entering stream tube illustrating two flow 


conditions into the nacelle. 
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in the internal flow. 

When Vi/Vo = 1, the flow is in an intermediate state 
between the other two conditions and no new limitations 
are imposed. This condition is of no interest, then, from 
the design viewpoint. 

Changes in local velocity — the geometric con- 
figuration of the nacelle, lip shape, bluntness of the 
forebody, etc., are perturbations on the flow conditions 
described above and are assumed to be minimized by the 
detailed design of the nacelle. No attempt is made here 
to define that detailed design except to assume that it is 
at least as good as the best of the empirical examples 
and data available to this investigation. 

A very frequently encountered set of requirements, 
as indicated by the literature, includes the following 
Situation: 

Cruise speed is approximately twice take-off speed. 
Cruise flow rate is approximately equal to take-off 
flow rate. 

Let the subscripts 1 апа 2 designate cruise and take-off 
conditions, respectively. If the inlet area is a constant, 

012Уі) Аі 7 Q5-Vio Ai 
Viz=Vi> 
but Voz=Vo2:2, so Vi/Vo]=%Vi/Vo2. 

If Vi/Vo,=0.7 to 0.8 as is E the case, 

(refs. 4,5,6,) then Vi/Vo2=1.4 to 1.6, but except when 


the inlet is very deeply submerged, these high velocity 
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ratios will cause cavitation inside the inlet, incurring 
high losses and perhaps choking the inlet. 

One way of decreasing this effect is to decrease the 
inlet velocity ratio at cruise by increasing the area. 
This approach very quickly reaches the minimum velocity 
ratio below which cavitation occurs in the external flow 
with its accompanying high drag and possible material 
damage. 

Às speed increases, the interval between these two 
limits becomes smaller, so the problem grows increasingly 
difficult as attempts are made to achieve higher speeds. 

A solution to the problem may be found in the use of 
. auxiliary inlets to increase the effective inlet area when 
high velocity ratios are required. 

If the existance of an auxiliary inlet is postulated, 
the velocity ratio at take-off is no longer coupled to 
the velocity ratio at cruise and both become design con- 
ditions as the total inlet area for each must be determined 
independently. 

An additional limitation is imposed by the elbow 
where the internal flow leaves the nacelle and enters the 
strut. This elbow will normally contain turning vanes to 
help minimize losses. Since the local velocities on the 
turning vane surface will be greater than the average 
velocity in the flow entering the elbow, a danger of 
cavitation exists. This turning vane cavitation is avoided 


by diffusing the flow in the nacelle to reduce the 
ge 





velocities and increase the static pressure. Therefore, 

a minimum diffusion ratio dictated by a maximum nacelle 
exit velocity, above which cavitation occurs in the elbow, 
is a third design limit. This limit must be checked at 
both cruise and take-off because it is a function of both 
the flow rate and the losses. 

In summary, then, both cruise and take-off are design 
conditions and cavitation avoidance establishes absolute 
limits on the minimum velocity ratio at cruise, the maximum 
velocity ratios at both cruise and take-off, and the 


minimum diffusion ratio. 
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Chapter 3. 
NACELLE SHAPE 

The principles of design used in jet engine nacelles 
and air intakes, already established in the field of 
aerodynamics, can be used for guidance in the selection 
of waterjet inlet nacelles (ref. 3 and 4). 

Küchemann and Weber (ref. 7) have developed the de- 
finition of a theoretical optimum nacelle contour by con- 
sidering the momentum balance between the free stream and 
the entry plane of a semi-infinite intake. The resultant 
force on the nacelle is given by 

Pie 3 ve A, (1 - V,/V,)?, 
‚but since there can be no resultant force on a semi-infinite 
body in potential flow, this force must be equal and 
opposite to the integral of the pressure over the forebody 


area of the nacelle. 


N An 
Ғы am) о, аА = -2/212 [1 - (V,/V,)* ] a4 
A; As 

If the pressure coefficient is a constant over the 
entire forebody then the local velocity, Va» is a constant 
and the projected frontal area, A,-A,, required to balance 
Fy is defined by the above equation. When the local 
velocity is limited to some maximum value to satisfy Mach 


number requirements or, in the case of waterjet nacelles, 


cavitation requirements, then the minimum projected frontal 
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area and Vas will satisfy the relationship 


AG tegen 


The optimum intake contour of Kuchemann and Weber is 
one which at its design Vi/V. achieves the goal of keeping 


Viax constant and at its specified value over the entire 


ma 
forebody and the equals sign applies. 

ET can be represented by the square root of its 
in the case of circular intakes. 
(1 - V;/Vo)2 1 2 
(Wmax/Vo)*-1 


It can be seen that a maximum value of D/Da is defined 


inverse, D, /D, , 


VD IN. 


-by a minimum value of Vi/Vo for any given V,. One of the 
design limits chosen in the previous chapter was avoidance 
of cavitation in the external flow at the relatively low 
velocity ratio at cruise speed. This can now be translated 
into a maximum value of D; /Dg. 

It is not practical to build such perfect intakes 
although intakes with reasonably uniform pressure dis- 
tributions at certain velocity ratios have been developed 
mer, 3, h, 7). In general, the velocity distribution 
will have peaks. 

The generation of practical contours with acceptiable 
pressure distributions has been largly a matter of 
experimental trial and error. Kuchemann and Weber de- 
Scribe three sets of geometrically similar shapes which 
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they designate as Classes A, B, and C. The geometric 
description of these are given in Table l and Figure 3. 

References 3 and A4 describe two sets of intakes 
whose pressure distributions were generated by numerical 
solutions to the potential flow problem. The Lockheed in- 
Bere chosen for input data for this study chiefly 
because the data covered a wide range of geometric variation 
and data on the interal flow distribution was available 
for each of the individual inlets. The major geometric 
parameters and peak pressure coefficients for each of the 
nine inlets is given in Tables 2 through ll. The external 
contour of these nacelles is very close to, but not quite, 
-elliptical. The internal lip contour is elliptical. The 
Kuchemann and Weber intakes are elliptical both internally 
and externally. Brown and Traksel state that this 
variation has an influence that is secondary to the length 
to diameter ratio, L/D , and D;/D,. Figure 4 shows the 
Similarity of the major parameters. The dashed line was 
tabulated as the basic geometric description of nacelle 
forebodies for SUBROUTINE NACEL. The three data points 
to the left represent intakes which had lower pressure 
coefficients than the others, so the dashed line represents 
the best of the available data points. 

To obtain a functional relationship for the variation 
of pressure coefficient with geometery, the peak pressure 


coefficient at each velocity ratio and at each angle of 
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attack was plotted against the value of DEDE and also 
against L/D,. L/D, was chosen as the appropriate parameter 
to use because it gave a smoother distribution of data 
points. The desired relationship was then obtained by 
drawing a smooth envelope that enclosed all points and 
passed through the points of highest C This envelope 
is considered to represent the "well designed" inlet 
mentioned in the problem statement based on the assumption 
that if a nacelle has already been designed which falls 
on that line (the discrete data points through which the 
line was drawn), then a nacelle could be designed for any 
other point on that line. The plots are shown in Figures 5 
through 10. Figures 11 through 16 show similar plots for 
the internal flow conditions. 

There can be no losses in the stream tube from out 
in the free stream up to the inlet plane, but from the inlet 
plane on into the ducting of the nacelle friction and 
turbulence take a toll of the head available in the internal 
flow. The amount of head loss is primarily a function of 
velocity ratio, Reynolds number and lip shape. Blackaby 
and Watson investigated experimentally several lip shapes 
including circular profiles, elliptical profiles and sharp 
profiles. Іп general, the elliptical profile caused less 
losses at a given Vi/V, than the others. The particular 
results for the lip which they designated 18E is given as 
a function of Reynolds number for several velocity ratios 


in reference 3, Figure 74 of Report No. 5. The loss 
=2 2 





coefficient, defined as ^ P/q., approaches a constant value 
for each velocity ratio for Reynolds numbers greater than 
8x10^. These constant values as read from the citation 
above are given in Table 12. The Ram Pressure Recovery, 


defined as 
RPR= (P = pa eo = р) - 1 = AP/q, 
is plotted in Figure 17, and points from that plot are 


tabulated in the subroutine. 
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Chapter 4 
THE DIFFUSER 

The inlet nacelle modeled in this study is assumed to 
consist of a streamlined body of revolution containing 
a nose intake opening, a diffuser and an auxiliary inlet. 
The forebody and intake lip serve to guide the flow into 
the internal ducting. That internal ducting includes the 
diffuser which has at least two primary functions. One 
is to slow the flow velocity to below the critical value 
for cavitation in the strut elbow. The other is to de- 
crease the flow velocity further to reduce ducting losses 
if it benefits the system as a whole. 
| Three main profile shapes of diffusers with circular 
cross section can be considered. They are bell, trumpet 
and straight walled conical diffusers. To avoid confusion 
of terms, a bell diffuser is one having its mean diameter 
greater than the mean diameter of a straight walled conical 
diffuser with the same inlet and exit areas while a trumpet 
diffuser has a mean diameter less than that of the cone. 
At similar inlet flow conditions, the bell shape results 
in a lower total head loss than the other two and the ` 
trumpet results in a higher loss (ref. 4, 9). 

In general, the diffuser with the least loss in total 
head is preferred for this application if other con- 
Siderations are equal because the difference in losses must 


eventually be made up by expenditure of prime mover power. 
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Current knowledge of boundary layer growth and be- 
havior in adverse pressure gradients does not allow accur- 
ate prediction of diffuser performance by analytic methods 
(ref. 3), and very little empirical data on bell diffuser 
performance that is directly applicable was found. Ref- 
erence 9 contains an experimental comparison of the three 
shapes which indicates that conical diffusers have pressure 
recovery characteristics very close to the bell and sig- 
nificantly better than the trumpet. Experimental results 
contained in references 9 through 13 can be applied to 
predict the losses of straight walled conical diffusers 
and can be considered a lower bound of that which can be 
expected from a bell. 

The total pressure loss for a conical diffuser can 
be predicted by 

Заа засе е 
where the second term is a ducting loss of a duct the same 
length as the diffuser and of a diameter equal to the mean 
diameter of the diffuser. For a cone, 
ат = (ру + 02)/2 . 
The first term accounts for diffusion losses as a percent- 
age of the losses in a sudden expansion. K, is the sudden 
expansion loss factor and is given by 
K sz (1- A/A))* . 


The proportionality factor C is a function of ex- 





pansion angle, 9, or the double angle, 20, and is determined 
experimentally as in reference 13. 

The minimum loss for any given diffusion ratio occurs 
when the value of 20 is in the vicinity of 59 to 109. Flow 
separation starts interfering with the flow and partially 
stalling the exit when 20 is in the vicinity of 209 
(ref. 11), therefore the region of interest is for 20 less 
than 209. In this region the curve of C vs. 20 can be 
approximated by 

C = 3.190 x 1072 (0)? + 8.452 x 107^ (6) 
where O is in degrees. 

There appears to be no reason in this application to 
decrease the angle below the value of minimum loss because 
the diffuser length increases and friction losses start to 
rise accordingly. For the computer model limits of 3° and 
10° were set on the minimum and maximum values of 8 which 
has the effect of establishing a maximum and minimum length. 

The desired total length of the nacelle is generally 
established by external hydrodynamic flow, structural and 
other considerations that do not include the diffuser and 
other internal length requirements. If the internal com- 
ponent lengths are less than the desired length, they will 
all fit within the EP If not, additional external 
drag due to the added length and wetted surface is generated 
and the power required to overcome that drag increment can 


be charged to the diffuser length since other internal 
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lengths are fixed. This becomes a criterion for determin- 


ing the.proper diffuser length. The sum of the power re- 
quired to overcome diffuser losses and the additional 
power due to the increment in drag is minimized with re- 
Spect to diffuser length. The calculations used in this 
determination are given in Appendix B. 

A limit on the minimum diffusion ratio is given by the 
requirement that the maximum diffuser exit velocity must 
be less than the critical velocity which will cause 
cavitation in the elbow. Experimental results in reference 
ll give exit velocity profiles for six diffuser examples, 
in two of which 26 = 10° and in four of which 20 = 12°. 
‘The diffuser area ratio was 2.0 in all cases, entry Mach 
number was 0.26, and boundary layer thickness ranged from 
0.003 to 0.03 times the entry radius. Integration of 
these distributions graphically gives a range of 


V - 1.15 to 1.50, with a mean value of 1.33 


max./ Vave. 

over the six cases. However, since boundary layer thick- 

ness has poor predictability, the average exit velocity is 

not known but it will be somewhat greater than the flow 

rate divided by the exit area. If Massa; is calculated by 
Vmax. = К Q/A5, 

k = 1.33 will probably give a low value and k = 1.50 may 

give a high value but will be on the safe side. 

The upper limit on V „x. will be the critical value 


of Vo determined from Appendix A. 
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Chapter 5. 
THE AUXILIARY INLET 

It was indicated in chapter 2 that an auxiliary inlet 
was probably required to obtain the flow rates necessary 
at take-off speed and remain subcavitating at cruise. 

The possible configurations for providing auxiliary 
inlets are many and varied. Slots in the nacelle wall with 
doors that open inward, leading edge openings on the strut 
and a movable nacelle lip are all possibilities, to mention 
only a few. 

The model for this study was assumed to be of the first 
type. Figure 3 shows the auxiliary inlet as dashed lines 
‘in the profile. It is assumed that the inlet is in the 
form of a series of Openings in an annular ‘ring around 
the nose opening with some fraction of that ring being 
structure to support the nose. The dimension X in the 
figure is the width of that annular ring in the direction 
normal to the inlet flow through the opening. 

The auxiliary inlet flow is assumed to enter and 
combine with the flow from the nose opening prior to enter- 
ing the diffuser because that is the area of lowest static 
pressure and, Since the inlet will only be used at take-off 
when the flow rate is high, it is important that the entire 
combined flow pass through the diffuser to prevent 
cavitation in the elbow. 

When the inlet is closed the length of ducting, Lys 


OD 





associated with it is considered to be a short length of 
pipe and the loss for that length of pipe is added to the 
other system losses. 

This type of inlet appears to be in much the same flow 
conditions as the “gliding scoop" described in reference 14. 
Since the auxiliary inlet is open only when the system de- 
mands a high flow rate its inlet velocity ratio is estimat- 
ed to be about 1.0 or slightly higher. In this condition 
the lip loss coefficient is approximately 0.2 (ref. 14). 
This value is the one used in the subroutine for auxiliary 


inlet losses. 
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Chapter 6. 
DESIGN RATIONALE 

The design of an inlet for the specific conditions 
imposed by the optimization program is based on comparison 
of the characteristics implied by the system parameter 
values with those stored in the data tables and sizing the 
individual components to remain cavitation free at the 
design conditions. The subroutine listing and a computer 
generated flow diagram are in Appendix C. This and the 
following chapter will attempt to clarify the reasoning 
used in its derivation. 

When the program calls SUBROUTINE NACEL, all physical 
constants, craft configuration, depth of submergence, 
free stream velocity, flow rate and angle of attack for 
each design point are known. Inlet velocity at cruise 
and D./Da are passed as system parameters., Actually, 
ШҰАР at cruise is the system parameter, but it is cal- 
culated from V; and а rather than passed. 

Initially, free stream static pressure, free stream 
dynamic pressure, incipient cavitation number, and total 


pressure for both cruise and take-off speeds are calculated. 


gu us 


a 2 
боі ` (ро 7 рул 


The length to diameter ratio, L/D , of the nacelle for 
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which the external minimum pressure coefficient is equal 
to the negative of the free stream incipient cavitation 
number at cruise is determined by interpolation in the 
data table where 


ya = (б, ИА a). 


The inlet to maximum diameter ratio, D;/D , correspond- 


m? 
ing to that L/D, is found in the tabulated geometric 

characteristics. The value obtained defines a nacelle for 
which the external flow on the surface is at the point of 


incipient cavitation and therefore iS a maximum value. 


(0/0) а = 200/0) 


If the Di/D, passed by the calling program is greater 


Ethan (D; /D ) that particular set of parameters is 


max? 
rejected by setting the nacelle weight to a very high 
number and returning to the calling program. If D/P, is 


less than (D,/D then it is cavitation free and the 


m) max: 
design proceeds. 
Inlet area is determined from the flow rate and inlet 
velocity at cruise which then allows Di and D. to be 
determined. 
The maximum allowed velocity ratio for both speeds is 
determined by interpolation in the data table again. 
ССА па 


This limit is imposed by cavitation on the inside of the lip. 


O., X) 


X E f(L/D,, Ol 


If at cruise, Vi/V, is greater than (Vi/V.) the 


max? 
parameters are rejected in the same manner as before. If 


com 





not, the take-off condition is checked to determine if 

an auxiliary inlet is required. The maximum flow rate 

allowed by (Vee s enters through the nose inlet, the 

remainder must enter through the auxiliary inlet. 

The pressure recovery of the lip is a tabulated 

function of V. /Vo, so the total pressure inside the inlet is 
iM er 

and the static pressure inside is given by 


2 
= m 
P; mí Pi e: NA e 
Pressure recovery of the auxiliary inlet is known 
from chapter 5 so the total pressure inside the auxiliary 
inlet is 
Faux = RPRaux o * Po ° 
The static pressure at a point is unique so the static 
pressure inside the auxiliary inlet is the same as that 
previously calculated from the flow inside the nose inlet. 
Therefore, the velocity through the auxiliary inlet is 
А 3 
же p [2 (Paus ü p,)/7] : 
The area of the auxiliary inlet is found by 
Aaux Е Qaux/ Vaux ° 
The total pressure of the. combined flow is estimated 


as the mass weighted average of the combining flows. 


P. = (Qaux Paux +.Q; P;)/ Q< 


where Qe = Qaux * Qi ° 


The average velocity of the combined flow is 
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calculated from the dynamic pressure 


Vi = Nep p p/p] * , 


and the pressure recovery coefficient to this point is 


E e Poa q. ° 
The diffuser entry and exit diameters are set by 
D, =D; , Do =0.9D,. 

The 0.9 factor is somewhat arbitrary. Some fraction 
of the nacelle diameter must be allowed for structural 
requirements and it was felt that ten per cent would be 
sufficient in most cases. 

The internal component length requirements are 
calculated in the manner of Appendix B and the sum of 
Bose lengths plus the lip length is used as the nacelle 
length subject to a minimum value of 5.5 - This factor 
is also arbitrary and is the average of the recommended 
values of Ly/D, found in the literature. Hoerner (ref. 15) 
shows from empirical data that the minimum drag nacelle 
on à wing has a L/D, value of approximately 3.0, but since 
the internal length requirements are always greater than 
3.0 (see Appendix B), it was felt that value might place 
an unncesssarily high drag penalty on the diffuser. 

The larger of the two flow rates, cruise or take-off, 
is used for sizing the diffuser. 

The diffuser loss, lip loss, and pipe loss, if any, 
are calculated for both cruise and take-off. Their sum in 


each case is used in the equation derived in Appendix A 
== 





to determine the critical velocity value at the diffuser 


exit for cavitation on the elbow turning vane. 


Ы 1 0ү2 1 # 5 
Verit p [65i tł- P1oss/3/2 V5 š Vf š Ov v) 
The maximum diffuser exit velocity as determined in 
chapter 4 is 
2max 102 Q/A,, : 


These two velocities are compared for both conditions. 


If e is greater than ae for either condition, the 
subroutine returns to the calling program as before with 
m arbitralily high weight value. If SNP is satisfactory, 
the design is complete. 

External drag is estimated by 

р = Cy À22V^ A 

where A is the external surface area of the nacelle. A is 
determined by dividing the nacelle in to three sections, 
the forebody from the nose to the point of tangency with 
a line parallel to the centerline at maximum diameter, a 
cylinderical midbody, and a tail section which is assumed 
to be symmetrical with the forebody. The forebody and 
tail can be approximated as paraboloids of revolution 


and their combined area is 


A, = (km /3)(D /2)2 [w - 1/(1 - W) 


p 
where M = [1 - h(2L/D )2] (ref. 4). 
For the cylinderical midbody, 
А Ер уо ли и 
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and then A = A + A. s 


The drag coefficient is from Hoerner (ref. 15) and 

is given by 
C, = Ce [1 + 1.5(0,/1,.)/? + 7(2, /142] ° 

The Schoenherr friction coefficient, Ce, is calculated 
by a function subroutine, CFS(R), in which the argument 
is Reynolds number, 

The weight estimate depends on materials and details 
of the structure which were not the concern of this pro- 
ject. Therefore, only a rough approximation is possible 
and the assumptions are certainly open to question. 

Internal ducting weight was specified to be propor- 
‘tional to the cross sectional area and the length such 
that W = 15.07 (Area) (Length) in pounds 
where area is in square feet and length in feet (ref. 1). 

The weight of the outer envelope is estimated by 
assuming that it is proportional to the surface area, 
the maximum wall thickness, and the specific weight. 

A maximum wall thickness of 0.05D, was used to be con- 
Sistant with the diffuser exit diameter assumption. 

The outer envelope of the nacelle, of course, is not 
solid but is made of inner and outer skins plus stiff- 
ening structure. The specific weight depends on the 
material which has not been specified. Two extremes of 
the weight spectrum which might be used for such a 


Structure are aluminum and steel. Some material with 
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a weight between these two is more likely and a material 
density of one half that of steel was used in the sub- 
routine. The actual weight in water may very well be 
close to zero if the voids are filled with oil or pitch. 

The proportionality constant, a measure of the vol- 
ume of the material in the wall compared to the volume 
of the wall itself, a packing factor, will be a small 
decimal on the order of 0.01 to 0.1. 0.03 was used in 
the subroutine. 

The resulting weight was increased by ten percent 
to allow for the auxiliary inlet mechanism. 

At this point the subroutine's job is completed and 
it returns to the calling program with the following 
quantities: 

a. the head loss summed from the free stream to the 

diffuser exit for both speeds, in feet, 

b. the drag of two nacelles for both speeds, in pounds, 

c. the exit area of two diffusers (one for each 

nacelle), in square feet, 

d. the weight of two nacelles, in pounds, 

e. the location of the vertical center of gravity of 

the nacelles below the keel of the craft, in feet, 

and 

f. the location of the longitudinal center of gravity 


of the nacelles forward of the transom, in feet. 
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Chapter 7. 
OFF-DESIGN PERFORMANCE 

Prediction of off-design performance may be carried 
out by SUBROUTINE NACEL only after a design has been com- 
pleted and returned to the calling program. The calling 
program must reset the value of a control variable named 
ISTRT from 1 which indicates design, to 3 indicating 
evaluation, and the value of a control variable named 
DELTA which has a finite value greater than 1l x 1077 
during design to O. which prevents rejection and return 
when cavitation limits are violated. 

The same initial values are known as when the design 
process started, and in addition, all values established 
during the design cycle are now stored for use during 
evaluation. 

The free stream conditions, total pressure, dynamic 
pressure and cavitation number, are calcuated for each 
Speed to be evaluated. 

The subroutine now proceeds through the same set of 
calculations as for design except that the auxiliary in- 
let, and diffuser sizing sections are omitted. Cavitation 
checks are made in the same manner, but the existing 
parameters are not allowed to be rejected. Instead of 
rejection and return a variable named CAV(I,J) has its 
value changed from O. to l. to indicate that cavitation 


is occuring at the evaluation speed indicated by I and at 
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the location indicated by J. J=l on the outside of the 
forebody, 2 inside the lip, and 3 in the elbow. 

On completion of an evaluation cycle the subroutine 
returns the value of the head loss in feet and the drag in 


pounds for two nacelles for each of the off-design speeds. 
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Chapter 8. 
CONCLUSIONS AND RECOMMENDATIONS 

Results of the subroutine's use with the optimiz- 
ation program are discussed in reference 1. In addition, 
to test it separately, the subroutine was called by a 
test program which stepped through values of D,/D, from 
О.д to 0.9 and values of cruise AE Proms TCO LOT 
three cruise speeds, 40, 45 and 50 knots. Take-off speed 
was assumed to be one half cruise speed. The required 
flow rate at cruise was arbitrarily set at 600 cubic 
feet per second with the flow rate at take-off ten per 
cent greater. The angle of attack was set at two degrees 
-and depth of submergence at five feet. 

Ram pressure recovery at the diffuser exit and cav- 
itation limits as calculated from the above input are 
displayed in Figures 18 through 21. The ram pressure 
recovery curves are unique to the particular flow rate 
combination assumed for this calculation and are present- 
ed as samples. The cavitation limits would be the same 
for any flow rate at the same angle of attack and depth 
of submergence. 

Several explicit and implicit assumptions were made 
in this study to simplify it or to try to bound a real- 
istic value. Further investigation might prove useful 
in any or all of the following areas: 


1. Very little information on bell diffusers that 
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could be directly applied in this study was found. 
Conical diffuser performance was assumed to be a lower 
bound on that expected from a bell. Information on bell 
diffusers similar to that in reference ll on other shapes 
would allow realistic prediction of their performance. 

2. The pressure coefficient information used from 
reference 3 was calculated rather than measured. Although 
the predicted values for simple cases (i.e., a sphere) 
was compared with experimental results, no mention is 
made of comparison of calculated and experimental data 
for a more complex shape such as an inlet nacelle. 

3. The drag calculation was made with the implicit 
assumption that the drag of a nacelle with a nose inlet 
and flow into that inlet is very close to the drag of 
that same body with no nose opening or inflow. Refer- 
ences 15, 16 and 17 support this assumption, but only a 
small number of cases with low D;/D, are examined. There 
is some indication in reference 15 Biss the range of 
ТАР for which this assumption holds decreases as D,/D,, 
increases. 

4. The effect of proximity to the free surface on 
drag was not included. For overall system performance 
comparison, the omission of this effect would probably 
not invalidate the comparison but there could be a slight 
bias in favor of the smaller nacelle in the case where 


there is a significant difference in size. 
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In its present form, SUBROUTINE NACELLE requires 
interpolation in data tables several times. Some of these 
tables could be approximated by functions to speed up the 
operation and reduce computer costs. The data for which 
this could be carried out include Table 12 and Figure 17. 
The data in Tables 5 through 16 presents a more difficult 
problem because the relationship is not single valued in 
all cases. Perhaps separating it into two or more regions 
with a different function applicable in each region offers 
a solution. Some unsuccessful attempts were made to ex- 
amine possibilities for the large data bank. 

One other simplification which would decrease com- 
‘puter time would be accomplished by elimination of the 
calculation of diffuser length. One of the results of 
the data run mentioned above was that the diffuser length 
calculated by iteration in the subroutine was nearly 
always approximately ninety per cent of the mean of the 
maximum and minimum allowed values. The equation 

Ly = 0.9(Lamin + Lamax)/2 


could be used to completely eliminate the iteration. 
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Table 1. Geometric Conficuration of Kuchemann and Webber 
Circular Intakes Designated Class A, B and C. 
* Dimensions of all three classes are defined by the 
following equations (adapted from ref. 7, p. 79). 

Do/Di = Ky 

Dm/2L = Kg - Kz (Ai/Am)3 

рі/211 - 1/K, (Ki - 1). 


Table l.a Values of the Constants. 


Class Кі К2 K4 К, 

А 1.15 0.2 12.5 15 

В 1.10 0.2 12.5 1.3 

C 1.03 0.2 [ло 1.0 
(ref. 7) 


Table 1.b Values of the Primary Ratios. 
Di/Dm L/Dm 


Class A and B Class C 
0.90 0.0732 0.1195 
0.80 0.1439 0.231 
0.70 0.309 0.452 
0.60 0.639 0.910 
0.50 1.263 1.575 
0.40 2.00 2.175 
0.20 2.1490 2.494 


*(See Fig. 3 for illustration of dimension symbols.) 





Table 2. Geometric Configuration of the Lockheed Intakes 


Nacelle *Do/Dm ш Di/Dm L/Dm 
A-10-550 0.40 9.0 0.333 2.20 
A-80-40 0.80 9.0 0.77 0.32 
A-50-280 0250 9.0 0.126 1.40 
A-75-60 0.75 9.0 0.713 0.45 
А-70-90 0.70 9.0 0.656 0.63 
A-55-280 0.55 9.0 0.483 lJ. 525 
A-60-135 0.60 9.0 0.541 O 288 
A-45-1400 0.45 9.0 0.369 1.80 
A-60-200 0.60 9.0 0.541 1.20 
(refe) 


*(See Fig. 3 for illustration of dimension symbols.) 
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Table 3. Variation of Minimum Pressure Coefficient with 
Velocity Ratio and Angle of Attack for Lockheed Nacelle 
A -40-550 

Internal Flow: 


Minimum Cp 


Vi/Vo e zQ 29 49 oe 

0.70 0.420 0.415 0.400 0.390 
0.80 0.260 0.250 0.240 0.220 
0.90 0.080 0.050 0.045 0.030 
1.05 -0.260 -0.255 -0.285 -0.550 
1.15 -0.500 -0.510 -0.550 -1.040 
1.25 -0.760 -0.765 -1.050 -1.630 


External Flow: 


0.70 -0.090 -0.110 -0.250 -0.440 
0.80 -0.090 -0.090 -0.160 -0.320 
0.90 -0.090 -0.090 -0.120 -0.230 
1.05 -0.090 -0.090 -0.115 -0.155 
2,15 -0.090 -0.090 -0.115 -0.145 
1225 -0.090 -0.090 -0.115 -0.140 
(ref. 3) 
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Table 4 


. Variation of Minimum Pressure Coefficient with 


Velocity Ratio and Angle of Attack for Lockheed Nacelle 


А-80-40 


Internal Flow: 


External Flow: 


Vi/Vo 


0.70 
0.80 
0.90 
1.05 
1.15 
1.25 


0.70 
0.80 
0.90 
1.05 
1.15 
1.25 


o. ZQ 
0.200 
-0.100 
-0.465 
-1.160 
-1.780 
-2.500 


-0.435 
-0 . 420 
-0.405 
-0.375 
-0.365 
-0.350 


Minimum Cp 


20 
0.050 
-0.285 
-0.700 
-1.425 
-2.220 


l -2.50 


-0.480 
-0.460 
-0 . 450 
-0.415 
-0 . 400 
-0.375 
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-0.540 
-0.500 
-0 . 490 
-0 . 440 
-0.425 
-0.410 
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-0.60 
-0.560 
-0.515 
-0.475 
-0 «460 
-0.440 


(ref.3) 





Table 5. Variation of Minimum Pressure Coefficient with 
Velocity Ratio and Angle of Attack for Lockheed Nacelle 
А-50-280 ° 


Internal Flow: 


=50= 


Minimum Cp 

Vi/Vo о 20 20 49 6? 
0.70 0.380 0.355 0.320 
0.80 0.200 0.160 0.120 
0.90 -0.015 -0.050 -0.100 
1.05 -0.370 -0.420 -0.485 
1.15 -0.640 -0.700 -0.770 
1.25 -0.940 -1.00 -1.080 

External Flow: 

0.70 -0.150 -0.160 -0.240 
0.80 -0.150 -0.155 -0.200 
0.90 -0.150 -0.155 -0.190 
1.05 -0.150 -0.155 -0.175 
1215 -0.150 -0.155 -0.175 
[25 -0.150 -0.155 -0.175 

(ref. 3) 





Table 6 


. Variation of Minimum Pressure Coefficient with 


Velocity Ratio and Angle of Attack for Lockheed Nacelle 


A-75-60 


Internal Flow: 


External Flow: 


Vi/Vo 


0.70 
0.80 
0.90 
1.05 
1.15 
25 


0.70 
0.80 
0.90 
1.05 
1915 
25 


a =0 
0.250 
-0.020 
-0.340 
-0.930 
01.440 
-2.050 


-0.360 
-0.350 
-0.340 


-0.320 


-0.310 
-0.300 


Minimum Cp 


20 
0.130 
-0.175 
-0.530 
-1.230 
-1.815 

-2.50 


-0.390 
-0.375 
-0.360 
-0.345 
-0.335 
-0.325 


ce 


1,0 
-0.020 
-0.350 
-0.780 
-1.600 
-2.370 


Z -2.50 


-0.445 
-0.410 
-0.395 
-0.370 
02565 
-0.345 
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Table 7. 


velocity Ratio and Angle of Attack for Lockheed Nacelle 


A-70-90. 


Internal Flow: 


Vi/Vo 
0.70 
0.80 
0.90 
1.05 
1.15 
25 

External Flow: 
0.70 
0.80 
0.90 
1.05 
1215 
1.25 


O =0 

0.305 
0.070 
-0.210 
-0.710 
-1.130 
-1.630 


-0.280 
-0.270 
-0.265 
-0.255 
-0.250 
-0.240 


Minimum Cp 


20 
O. 215 
-0.040 
-0.350 
-0.930 
-1.420 
-2.120 


-0.310 
-0.295 
-0.285 
-0.2'70 
-0.265 
-0.260 


-52- 


1,0 

0.100 
-0.180 
-0.535 
-1.225 
-1.920 


2 -2.50 


=О > 
-0.335 
-0.315 
-0.300 
-0.290 
-0.280 


Variation of Minimum Pressure Ое Оол vith 
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-0.455 
-0.405 
-0.365 
-0.330 
-0.320 
-0.310 


(ref. 3) 





Table 8. 


Velocity Ratio and Angle of Attack for Lockheed Nacelle 


A-55-280. 


Internal Flow: 


Vi/Vo 
0.70 
0.80 
0.90 
1.05 
1.15 
1.25 

External Flow: 
0.70 
0.80 
0.90 
1.05 
1.15 
125 


e zQ 
0.415 
22286 
0.030 
-0.330 
-0.595 
-0.880 


-0.125 
-0.125 
-0.125 
-0.125 
-0.125 
-0.125 


Minimum 


20 
(95255 
0.200 

-0.010 
TOPENG 
-0.650 
-0.950 


-0.225 
-0.145 
-0.135 
-0.130 
-0.130 
-0.130 


57 


Cp 
1,0 
0.350 
0.160 
-0.070 
-0.440 
-0.780 
-1.380 


-0.450 
-0.270 
-0.185 
-0.170 
-0.150 
-0.150 


Variation of Minimum Pressure Coefficient with 
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-0.700 
-0.500 
-0.320 
-0.205 
-0.190 
-0.170 


(ref. 3) 





Table 9. 


Velocity Ratio and Angle of Attack for Lockheed Nacelle 


A-60-135 


Internal Flow: 


Vi/Vo 
0.70 
0.80 
0.90 
1.05 
15 
25 

External Flow: 
0.70 
0.80 
0.90 
1.05 
155 
1925 


ce zQ 

0.300 
0.080 
-0.175 
-0.620 
-0.965 
-1.340 


-0.255 
-0.255 
-0.250 
-0.245 
-0.240 
-0.240 


Minimum Cp 


20 
9227900 
-0.010 
-0.280 
-0.755 
-1.120 
O 


-0.280 
-0.275 
-0.270 
-0.265 
-0.265 
-0.260 


a 


1,0 

0.140 
-0.115 
-0.405 
-0.910 
-1.330 


£-1.80 


-0.315 
-0.300 
-0.295 
-0.280 
-0.280 
-0.275 


Variation of Minimum Pressure Coefficient with 


60 
0.040 
-0.235 
-0.560 
-1.135 
-1.620 
-2.230 


-0.380 
-0.350 
-0.330 
-0.310 
-0.300 
-0.295 


(ref. 3) 





Table 10. 


Variation of Minimum Pressure Coefficient with 


Velocity Ratio and Angle of Attack for Lockheed Nacelle 


A-45-400 


Internal Flow; 


Vi/Vo 
0.70 
0.80 
0.90 
1.05 
5 
25 

External Flow: 
0.70 
0.80 
0.90 
1.05 
15 
1.25 


© =0 
0.400 
0.230 
0.025 
-0.320 
-0.570 
-0.850 


-0.120 
-0.120 
-0.120 
-0.120 
-0.120 
-0.120 


Minimum 


20 
0.385 
0.205 
0.000 

-0.340 
-0.600 
-0.880 


-0.140 
-0.130 
-0.130 
-0.130 
-0.130 
-0.130 


25568 


Ср 
4,9 
0.370 
0.185 
-0.020 
-0.365 
-0.630 
-0.915 


-0.310 
-0.210 
-0.155 
-0.150 
-0.150 
-0.150 


60 


(ref. 3) 





Table 11, 


Velocity Ratio and Angle of Attack for Lockheed Nacelle 


А-60-200 


Internal Flow: 


Vi/Vo 
0.70 
0.80 
0.90 
1.05 
19.15 
1.25 

External Flow: 
0.70 
0.80 
0.90 
1.05 
1.15 
N25 


O =O 

0.400 

0.205 
-0.015 
-0 . 400 
-0.685 
-1.000 


-0.160 
-0.160 
-0.160 
-0.160 
-0.160 
-0.160 


Minimum Cp 


20 
0.360 
0.150 

-0.080 
-0.485 
-0.800 
-1.265 


-0.200 
-0.175 
-0.160 
-0.160 
-0.160 
-0.160 


"56 


1,0 
0.300 
0.075 

-0.165 
-0.620 
-1.140 
-1.890 


-0.385 
-0.240 
-0.210 
-0.185 
-0.180 
-0.180 


Variation of Minimum Pressure Coefficient with 


6° 


-0.660 
-0 . 420 
-0.280 
-0.230 
-0.215 
-0.210 


(ref. 3) 





Table 12. Pressure Loss Coefficient for an Elliptical 


Profile Lip 


Vi/Vo АР/до 
0.6 (215 
0.8 0.033 
1.0 0.043 
1.2 0.085 
leh 0.130 
1.6 0.173 


Reynold's number=8 X 104 


(Adapted from Fig. 74, Report No. 5, ref. 3) 


p = 
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Figure 4. L/Dm vs. Di/Dm for the Lockheed intakes and 


Kuchemann and Weber intakes Classes A, B and C. 
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Figure 12. Minimum Cp vs. L/Dm in the internal flow at 
Vi/Vo - .80 
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Figure 13. Minimum Cp vs. L/Dm in the internal flow at 


Vi/Vo - .90 
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Figure 15. Minimum Cp vs. L/Dm in the internal flow at 


Vi/Vo = 1.15 
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Figure 17. Ram Pressure Recovery vs. ШАР for the elliptic- 


al profile lip (Blackaby and Watson 18E). 
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Figure 18. Ram Pressure Hecovery at the diffuser exit vs. 


ШАЙ for a 40 knot cruise speed 
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Figure 19. Ram Pressure Recovery at the diffuser exit vs. 


V/V. for a 45 knot cruise speed 
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Figure 20. Ram Pressure Recovery at the diffuser exit vs. 


АЙ for a 50 knot cruise speed 
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Figure 21. Cavitation limits as a function of D./D vs. 
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Appendix A. CAVITATION CRITERA 
The total pressure head at any point in the system can 
be determined with respect to the free stream head by trac- 
ing along a stream line from out in the free stream to 
the point of interest. 
H = Ho - Hioss 
op P P. pd oscmNESbmS of pressure: 

Using Bernoulli's equation and allowing the static 
head term to reflect all changes in elevation or 
submergence, 

5/28 =h- V*/2g s 
ог р + 2) vé = р + ¿PV + ав" 
The cavitation number, gp, сап be defined as 
C 8 (Pret. = p)/5⁄ V“ 


and VET are the reference pressure and 


ref: 
Where D. 


ef. 
velocity to which the cavitation number is related and 
p is the local static pressure at the point of interest. 
If the free stream conditions are used for reference, 
then pg апа У. аге the reference pressure and velocity. 
In this case 
a = 1 2 _ 1 2 
Po TE zi 2 15 + Ploss 
so 07 = (р. - р) /2 212 = V2/V2 - 1 + P. 42/3 2х2 
0 Po ` P//2/ 15 о loss/?^ *g 
Cavitation is assumed to occur when the local static 
pressure is equal to the vapor pressure of water. The 
incipient free stream cavitation number is defined by 


por Ung py) /3^ V6 Я Td ү5)-1 ii P1oss/ 37 VS. 
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This definition allows the prediction of cavitation using 
the free stream incipient cavitation number if the losses 
are known. 

Consider the situation in Figure (3). The diffuser 
entry is designated station l, and the diffuser exit 
station 2, to remain consistant with the rest of the text. 
The turning vane is not intended to be a particular type 
or design, but is expected to be well designed and sized 
for the elbow and diffuser with which it is used. The 
sketch is conceptual, only. The subscript tv will designate 
conditions on the turning vane. 

To determine the cavitation number on the turning vane, 
conditions at station 2 can be used and, in turn, related 
to the free stream. As before, 


e 


EV Р 11055 


2—tv 

If the turning vane starts near station 2, then any 
loss between there and the point of lowest pressure can 
be neglected as it will be very close to the leading edge. 

Pey = P2 Pey + 2P Vey = Po * Ир 
The cavitation number on the turning vane is, 
a 1 2 

When Ру АЕ cavitation starts, so the incipient 


cavitation number is 
Е 2 
Оъуі = (рә - ру) /2 215: 


Ő tyi is a characteristic number determined from 
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experiments and is therefore known (ref.28). The con- 
ditions at station 2 which will cause cavitation on the 
turning vane can be determined. Solving the equation 
above for Pos 

Po = бұ, à ovs t Pys 
but P can be ЕА in terms of the free stream 


conditions, which are known. 


2 2 
> mE - 
Ро * a 7 eee 


P5 E: 
2 2 2 
a = 1 _ à = 
or бу 8 /^ V, tape Dot UE г РУ, CE 
P in this expression is the sum of all losses from the 


loss 


free stream to station 2. V, is the only variable re- 


maining which indicates the conditions at station 2, so 


solving for Е. 


2 
V (2 + с 


= = 1 2 E 1 2 2 
tvi! | (5, ру)/2 2 10 БЕ P. oss 57 V. ] % 


4 


Тһе first term in the brackets is the free stream 


incipient cavitation number. The value of V, which gives 


2 
rise to incipient cavitation on the turning vane is 


2 1 1 
= ^ 1 2 2 
Vo Ж + 1 Pioss/? M | v (1 + ШО) ° 
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Appendix B. GEOMETRIC DERIVATIONS 

There are three variable lengths in the nacelle which 
are not a function of the forebody geometry and which must 
be calculated in each design. They are the diffuser length, 
Lis the additional length due to the auxiliary inlet, Ly, 
and a length after the diffuser exit containing part of 
the elbow which directs the flow upward into the strut, the 
strut inlet length. 

1. The Strut Inlet Length 

Assume that the thickness to chord ratio of the strut 
all along its length from the nacelle to the free surface 
is determined by the criterion 

IC - q + — 1] PELIS (ref. 1,4) 

where Oi is the local free stream incipient cavitation 
number. Also assume an elliptical cross section for the 
strut and a rectangular cross section for the internal 
ducting (ref. 1). 

The respective areas are 

Astrut = (1/4) tm o 
Aduct = 5 tm Co (reta EE 

Aduct must be the same as the exit area of the diffuser, 
А2. 

To allow for structure апа fairing, assume that the 
desired additional nacelle length is close to the value of 
the strut chord, со. 


242 = 0.12 c,*, so that 


For Os = 0295 tm - 0.12 Co» O 3 


ES 





Co = 3.515000 E 

2. Length Due to the Auxiliary Inlet 

If the auxiliary inlet has the configuration shown 
in Figure 3, the inlet area is a section of a right circ- 
ular cone with a minimum diameter р;, and slant height X. 

The angle of inclination from the centerline of that 
cone section, shown in the figure as 6,, is assumed to be 
formed between the centerline and a line passing through the 
nose leading edge to the point of tangency of the forebody 


Ф 


at maximum diameter, D, 


e, = tan-l [(D, - D;)/2L] 

The mean diameter of the conical section is then 

D = D; + X sin 9, 
and its area is 
ТА = TDX . 

The area of the inlet will actually be some percentage 
of the area above because some allowance must be made for 
structure. 

sa - KrD X 
K = 0.8 is used in the subroutine. 


Finally, solving for X in the above expression, 


"d 
X = (0.5/sin 04) [(D;* + 1.273 An y Sin 0,/K)2- Di] , 


and D X COSMOS me 
3. Diffuser Length 
The length of the diffuser, Lg» becomes a factor in 


the determination of nacelle length when the sum of the 
2р р 





internal lengths exceeds the minimum length, 2.50 Тһе 
excess length is chargeable to the diffuser because all the 
other internal length requirements are fixed by other con- 
Siderations. The difference in length is 
Ly = Liip + L. + La + 3.5k5 D> - 5.5007 - 
The change in drag is given by 
A Drag = Cp iv^ INC 

because Cy is essentially a constant for small changes in 
Ly + 

The change in length can be considered to be all in a 
cylinderical section of diameter D. 

AS MD s 
The increase in power due to the increase in drag is 
Lower - 0) 1 2 Vo? Tr Dg Alyn € 
The power to overcome the diffuser losses is found by 
A Powers = (CK, « fLlg/d4)À2V;^Q . 

The total increase in power required to operate the 

diffuser is then 
A Power = APowerz + A Powerz . 
In the subroutine, the minimum value of A Power is 


determined by iteration on lg. 
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Appendix C. SUBROUTINE LISTING 

SUBROUTINE NACEL and its flow diagram listed in the 
following pages are in the form for which the rationale is 
described in chapters 6 and 7. Similar listings for 
FUNCTION TABLE, which was not programmed by the author but 
was furnished for common ues with reference 1, follow 
NACEL. FUNCTION CFS is also shown for completeness but no 
flow diagram was generated for it because of its simplicity. 

The FORTRAN names below, used for various variables and 
constants in the subroutine, are described as an aid in 
understanding the logic and specific calculations. The 
symbols are listed in approximate order of first occurance 
beginning with the first executable statement. 

ZK - A constant decimal less than one which in- 

dicates how much of the annulus occupied by 


the auxiliary inlet is available for opening. 


SPO - Static pressure in the free stream in p.s.f. 

HS - Depth of submergence in feet. 

HA - Pressure head of the atmosphere in feet of 
water. 

RHOW - Density of water in slugs per cubic foot. 

G - Acceleration of gravity, 32.174 feet/sec“ 

PVP - Vapor pressure of water in p.s.f. 

PV - Vapor pressure of water in feet of water. 


VELR(1) - Inlet velocity ratio at speed I. I =1 


indicates cruise. I = 2 indicates take-off. 
=o = 





VI(I) - 
VO(I) 
SIGTV - 


JNUMB - 
ISTRT - 


NUMB - 


QO(I) - 


SIGI(I) 
PTO(I) - 


I greater than 2 indicates a performance 
evaluation speed. 

Inlet velocity in feet/sec. 

Free stream velocity in feet/sec. 
Incipient cavitation number on the turning 
vanes referenced to diffuser exit pressure 
and velocity. 

À counter. 

A control variable indicating whether the 
current calculation is part of the design 
process (ISTRT = 1), or part of the eval- 
uation process (ISTRT = 3). 

A control variable indicating the total 
number of speeds to be considered. 

Free stream dynamic pressure at speed I in 


реБ 


- Free stream incipient cavitation number. 


Free stream stagnation pressure in p.s.f. 


TRIM(I) - Angle of attack in degrees. 


CPEX - 
CDUMX - 
CDUMY - 


External minimum pressure coefficient. 
A dummy array generated during interpolation. 


A dummy array generated during interpolation. 


CPEXT(J,K,I) - Tabulated minimum pressure coefficients 


in the external flow at the Jth velocity 
ratio, Kth length to diameter ratio, and Ith 


angle of attack. 
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CPINT(J,K,I) - Tabulated minimum pressure coefficients 
in the internal flow. 

VRT - Tabulated array of velocity ratios. 

IL, KL, ML, NL - Control variables indicating the 
order of interpolation desired. 

XDT - Tabulated array of angle of attack. 

XDTT(I) - A dummy array generated during interpolation. 

DIDMT - Tabulated array of Di/Dy. 

DIDMX - The maximum value of D,/D, for a V;/V, as 
determined from interpolation in the data 
table. 

XD - A working variable for length to diameter 
ratio. 


WGTS(1,1) - Wet weight of the nacelle in pounds. 


DIDM - Current value of Di/Dg. 

QI - Flow rate at cruise through one inlet in 
с.г.5. 

AI - Nose inlet area in feet squared. 

DI - Current value of D; in feet. 

DM - Current value of D, in feet. 

ELEXT - Length of the forebody in feet. 


VRTT(I) - A dummy velocity ratio array generated during 
interpolation. 

VRTEX(I) - A dummy velocity ratio array generated during 
interpolation. 


VRMAX(I) - Maximum allowed velocity ratio for a given 
“аса 





D;/Dg as determined from the tables. 
QIN - Working variable for flow rate through the 


nose inlet at take-off in c.f.s. 


QC - Required flow rate at take-off in c.f.s. 

QAUX - Flow rate through the auxiliary inlet in 
OILS 

KDEX - A counter. 

VI2 - Working variable for VI(2). 

PRLT - Tabulated array of pressure recovery coeffic- 


ients for the lip. 


PRLT2 - Lip pressure recovery coefficient at take- 
OTT 

PTI - otagnation pressure in the internal flow in 
p.s.f. 

SPI - Static pressure in the internal flow in p.s.f. 

AIAUX - Area of the auxiliary inlet in feet squared. 

VIAUX - Auxiliary inlet velocity ratio. 

PTAUX - Stagnation pressure inside the auxiliary 
inlet in p.s.f. 

PRAUX - Pressure recovery coefficient of the aux- 
iliary inlet. 

DYP - Dynamic pressure inside the auxiliary inlet 
in p.S.i. 

PC - otagnation pressure in the combined flows 


from the auxiliary and nose inlets in p.s.f. 


DLIP(I) - Inlet losses at speed I in p.s.f. 
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QDIF - Flow rate through the diffuser iN Cl oL 


ELENT - Interior length of the lip in feet. 

PHI - The angle labeled O, in Figure 3. 

ELAUX - Interior length due to the auxiliary inlet 
in feet. 

ELMAX - Maximum permissible length of the diffuser 
in feet. 

ELMIN - Minimum permissible length of the diffuser 
in feet. 

QDIF - Flow rate through the diffuser inc.f.s. 

EL - Current value of the diffuser length in feet. 

DEL - Change in diffuser length in feet. 

ELN - Length of the nacelle in feet. 


ELD, ELL - Working variables for ELN. 


DDM - Mean diameter of the diffuser in feet. 
XKT - Sudden expansion loss coefficient, K}. 
REL - Reynolds number based on length. 

RED - Reynolds number based on diameter. 

DM - Ratio, Dm/ Ly. 

CDRG - Drag coefficient. 

ANGL - Deffuser half angle of expansion, labeled 


O in Figure 3. 

POW - Current value of the power charged to the 
diffuser in ft.lbs./sec. 

POWI - Previous value of POW. 

ELDIF- Final value of diffuser length in feet. 
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REND - oame as RED. 


RENL - oame as REL. 

DDIF - Diffuser loss in p.s.i. 

PLOSS - oum of diffuser, pipe and lip losses in p.s.f. 
VAOUT - Average diffuser exit velocity in ft./sec. 
VCRIT - Critical diffuser exit velocity at which 


cavitation occurs on the turning vanes in 
IE cece 

VMAX - Maximum diffuser exit velocity in ft./sec. 

DELH(I,1) - PLOSS converted to head in feet of water 
at speed I. 

PRC(I) - Ram pressure recovery coefficient at the dif- 
fuser exit. 

CD(I) - Drag coefficient. 

AEXN - Wetted area of the nacelle in feet squared. 

CPOD(I) - Drag of two nacelles in pounds. 

AREA(1) - Diffuser exit area of two diffusers in feet 
squared. 

CGS(1,1) - Distance of the vertical center of gravity 
of the nacelles from the keel in feet. 

CGS(2,1) - Distance of the longitudinal center of 
gravity of the nacelles from the transom in 
feet. 

CGS(3,1) - Distance of the VCG of the contained water 
from the keel in feet. | 


CGS(4,1) - Distance of the LCG of the contained water 
-88- 





from the transom in feet. 


RHOD - Density of steel in slugs per cubic foot. 
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Appendix D.  ALTERNATE DESIGN RATIONALE USING ONLY VELOCITY 
RATIO AS A SYSTEM PARAMETER 

Another,approach to the design of the inlet nacelle 
can be made using the maximum allowable Di/D. at a given 
ТАА in stead of allowing it to vary as a Parane eroe 
design proceeds as previously described except that the 
diameter ratio determined from interpolation in the data 
table is used to calculate the nacelle geometry. If, after 
the diffuser is sized, the exit velocity is too high, then 
instead of rejecting that set of parameters, a new value 
of рр based on the maximum allowed exit velocity is cal- 
culated and the subroutine goes back to the beginning of 
the design phase where a new design begins with this value. 

Performance evaluation proceeds exactly as before. 

The following pages contain a listing of the subroutine 


in this form with a flow diagram. 
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Appendix E.  ASSUMED CONSTANTS AND LIMITS 

| lt was necessary to make numerous assumptions and 
impose limits on various quantities in the subroutine. Some 
of the assumptions were somewhat arbitrary and may be mod- 
ified whenever good reasons to do so are discovered. 

Some of those assumptions and limits are discussed 
here in the approximate order in which they are encountered 
in SUBROUTINE NACEL. 

The first assumption that should be mentioned 
is that everything in the model is based on a nacelle with 
a circular cross section. No attempt was made to account 
for other shapes. 

The decimal constant, ZK, which determines the 
amount of available inlet area in the auxiliary inlet ring 
Баз ап arbitrary value of 0.8. The real case may vary 
widely from this depending on the details of the design. 

The incipient cavitation number on the turning 
vanes, SIGTV, has a value of 0.4. This value was accepted 
as a reasonable one for the speed range of interest. No 
examples with which to compare were found in the literature. 

The angle of attack at cruise, TRIM(1), has been 
artificially limited to a maximum of three degrees because 
the external pressure coefficients are notsingle valued 
with respect to L/D, or D¿/D, for a> 4°. The highest 
tabulated value for angle of attack is 6° and extrapolation 


would probably not be valid. 
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The range of tabulated velocity ratios is from 
ӘЙ со 1.25: Extrapolation is allowed in the subroutine 
but the results could contain large errors if extrapolated 
ux far. 

The largest value of L/D, in the tables is 2.50. 
This is also allowed to be extrapolated and the same remarks 
apply as for the velocity ratios. 

D. /D, is limited to 0.90, which assumes a wall 
thickness allowance of 0.05D,. This value is arbitrary and 
may be modified as desired. 

The pressure loss factor for the auxiliary inlet 
is assumed to be 0.2 which contains an underlying assumption 
that its inlet velocity ratio is always near 1.0. 

Diffuser lengths are limited to those correspond- 
ing to double angles of 6% to 20%. See chapter 4 for an 
explanation of the limits. 

In the wetted surface calculation, the area of 
fairings into the foils and strut, as well as the missing 
area in way of the foils and strut are ignored. The nacelle 
is assumed to be a solid body of revolution. 

Interference drag between the nacelle and the 
other components is assumed to remain constant even though 
the nacelle size changes. No attempt is made to account 
for the effect of V./V. or the depth of submergence on drag. 


Chapter 8 contains recommendations concerning this subject. 


Ze 





The weight estimate, as discussed in chapter 6, 
is only a rough approximation with a series of assumptions 


as to the construction and material. 


-134- 














Thesis 
C713 


Conner 
A study of ram tvpe 


intake parameters, 





VLLL 
3 2768 002 09328 8 


M 
DUDLEY KNOX LIBRARY 






































